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We define and study a variant of the center of mass of a poly- 
gon and, more generally, of a simplicial polytope which we call the 
^ Circumcenter of Mass ( CCM). The Circumcenter of Mass is an affine 

(-H combination of the circumcenters of the simplices in a triangulation 

of a polytope, weighted by their volumes. For an inscribed polytope, 
CCM coincides with the circumcenter. 
i_i Our motivation comes from the study of completely integrable dis- 

crete dynamical systems, where the CCM is an invariant of the discrete 
j> bicycle (Darboux) transformation and of recuttings of polygons. 

We show that the CCM satisfies an analog of Archimedes' Lemma, 
a familiar property of the center of mass. We define and study a gen- 
O eralized Euler line associated to any simplicial polytope, extending 

the previously studied Euler line associated to the quadrilateral. We 
<^ show that the generalized Euler line for polygons consists of all cen- 

ters satisfying natural continuity and homogeneity assumptions and 
• • Archimedes' Lemma. 

. !_h Finally, we show that CCM can also be defined in the spherical 

^ and hyperbolic settings. 

1 Introduction 



This note concerns the following geometric construction. Let P = (Vi, V2, . . . , 
V n ) be an oriented n-gon in the plane and let O be a point not on the sides of 
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P (or their extensions). One has a triangulation of P by n oriented triangles 
OViVi + i where the index is understood cyclically mod n; let Cj be the circum- 
center of the ith triangle and Ai its signed area. Then the sum A\ + • • • + A n 
is the signed area A(P) of the polygon; we assume that A(P) ^ 0. Consider 
the weighted sum of the circumcenters 

n . 

CCM(P) = J2A(k) Q - (1) 

We call this point the Circumcenter of Mass of the polygon; we shall show 
that it does not depend on the choice of the point O. 




Figure 1: Construction of the Circumcenter of Mass 



In other words, we break P into triangles OV^+i and replace each such 
triangle with a point-mass equal to its area and located at its circumcenter. 
Then CCM(P) is the center of mass of these point-masses, see Figure [l] 

We learned about this construction from [1] where this center is mentioned 
as a point invariant under the transformation of the polygon called recutting. 
In [T7] , we proved that the Circumcenter of Mass is invariant under another 
map, called the discrete bicycle (or discrete Darboux) transformation; see 
[101 ITS] . Although our motivation comes from the study of these discrete 
dynamical systems, we believe that the Circumcenter of Mass is an interesting 
geometric construction and it deserves attention on its own right. 

If one replaces the circumcenter of a triangle by its centroid then, instead 
of CCM(P), one obtains the center of mass of a polygon P (more precisely, 
of the homogeneous lamina bounded by P), which we denote by CM(P). 
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In the case of CM(P), one may safely allow degenerate triangles, whose 
centroids are well defined and whose areas vanish, and which consequently do 
not contribute to the sum. However, not all degenerate triangles are allowed 
in the definition of the Circumcenter of Mass: although their areas vanish, 
their circumcenters may lie at infinity, in which case formula ([I]) fails to make 
sense. This discontinuity with respect to the location of point O is a subtle 
point of the definition of CCM(P). 

The contents of this paper are as follows. In Section [2j we prove that the 
definition of the Circumcenter of Mass is correct, that is, does not depend 
on the choice of the point O. 

In Section [3j we establish a number of properties of the Circumcenter of 
Mass: 

1. The circumcenter of mass satisfies the Archimedes Lemma: if a poly- 
gon is a union of two smaller polygons, then the circumcenter of mass 
of the compound polygon is a weighted sum of the circumcenters of 
mass of the two smaller polygons. As a consequence, one may use any 
triangulation of P to define CCM(P). 

2. If P is an equilateral polygon then CCM(P) = CM(P). 

3. One can naturally extend the notion of the Circumcenter of Mass to 
smooth curves. However, this continuous limit of the Circumcenter of 
Mass coincides with the center of mass of the lamina bounded by the 
curve. 

We start Section [4] with the observation that one can replace the circum- 
center or the centroid of a triangle by their affine combination and define 
a new "center" of a polygon. The affine combinations of the circumcenter 
and the centroid of a triangle form the Euler line of the triangle. This line 
contains a number of other interesting points, for example, the orthocenter. 
Applying our triangulation construction to a polygon P yields a line that we 
call the generalized Euler line of P. 

We investigate whether there are other natural assignments of "centers" 
to polygons, satisfying natural properties: the center depends analytically on 
the polygon, commutes with dilations, and satisfies the Archimedes Lemma. 
We show that all such centers are points of the generalized Euler line. 

We also show that the Euler line is sensitive to symmetries, both Eu- 
clidean (reflections, rotations), and others, such as equilateralness. 
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In Section |5j we show that the construction of the Circumcenter of Mass 
extends to simplicial polyhedra in IR n in a rather straightforward manner. 
In Section [6j we extend our constructions to the spherical and hyperbolic 
geometries. 

Whenever possible, we try to provide both geometric and algebraic proofs, 
so some statements in this paper are given multiple proofs. 

A historical note. After this paper was finished, we learned from B. 
Griinbaum that the notion of the circumcenter of mass appeared in the book 
by Laisant [11] under the name "pseudo-centre"; the discovery is credited 
to the Italian algebraic geometer Giusto Bellavitis. B. Griinbaum and G. 
C. Shephard obtained some of the result that appear below in 1993, but 
they did not publish this material. We are grateful to B. Griinbaum for this 
information. 



2 Correctness of definition. Degenerate tri- 
angles 

Let us prove that the Circumcenter of Mass is well defined. 

We shall use the following notation for the CCM of two polygons. If P 
and Q are polygons, define 

It is easy to check that this operation is commutative and associative. We 
also define 

rrui p « n\ - A{P) CCM ( p ) Z A (Q) CCM(Q) 
CCM(PQQ)- A{p) _ A{Q) . 

We consider CCM(QP) to have the same CCM as P but with mass of 
opposite sign. 

Lemma 2.1 Given a quadrilateral ABCD, one has: 

CCM {ABC ® CD A) = CCM (BCD © DAB). 
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Proof. Let Ca, C b , C c and C D be the circumcenters of ABCD, AACD, 
etc., see Figure [2| It suffices to see that CaCcC\CbCd, divides the diagonals 
of CaCbCcCd at the ratio corresponding to the areas of the triangles. We 
show this for one pair, and by symmetry, it will follow for the other. 




Figure 2: Proof of Lemma 2.1 



An easy angle count shows that ZCbCdCd = ZBAC, ZCaCdCb = 
ZACB, ZCcCbC d = ZCAD, ZC A C C C D = ZABD, ZC C C A C D = ZCBD. 
Indeed, the angles CaCdCb and ACB have pairwise perpendicular sides, see 
Figure |2j and likewise with the other pairs of angles. 

To see that the CCM divides diagonal C A C C in ratio A(ABD) : A(CBD), 
it is enough to show that A(ABD) : A(CBD) = A(C D C C C B ) : A(C D C A C B ). 
Using the angles from above and the Law of Sines, 

A(C D C C C B ) _ C C C D ■ C B C D sinjZBAC) _ C c C D sin{ZBAC) _ C C C D ■ BC 
A(C D C A C B ) ~ C A C D ■ C B C D sin(ZACB) ~ C A C D sin(ZACB) ~ C A C D ■ AB 
BCsin(ZCBD) BC ■ BDsin(ZCBD) A(ABD) 



AB sm(ZABD) AB ■ BD sin(ZABD) A(CBD) ' 



as needed. □ 



We use Lemma 2.1|to prove that the Circumcenter of Mass is well defined. 



We wish to show that CCM(©" =0 1 ^l / j + iO) is independent of O. We proceed 
by induction on n, the number of vertices of the polygon P. 
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The base case is n = 3. By Lemma 2.1 



ccMiy^o © v x v 2 o © v 2 v o) = ccMiy^o © v 2 v x o © ov v 2 ) 

= CCM(V Q VxO © ViOF © V V 2 V l ) = CCM(eV V 2 V 1 ) = CCM{y^V x V 2 ). 

Now assume the inductive hypothesis for n — 1 > 3 and consider the case 
for n: 

CCM(®toViV i+ iO) = CCMd^ViVi+iO) © V n ^ 2 V n ^O © V n ^V O) 



CCM({®^ViV i+l O) © K-1K-2O © OVbK 



n— 1 , 



By Lemma [2. 1[ this expression is equal to 

CCM{{@^ViV i+1 0) © K- 2 Oy © HK-1K-2) 
= CCM((®^ViVi +1 0) © K-2W) © CCM(V n - 2 V n ^V ). 

The first term of the last expression is the CCM of a closed (n — l)-gon. By 
induction, it is independent of O. The second term is clearly independent of 
O. This completes the induction. 

Next, we deduce coordinate expressions for the Circumcenter of Mass. 
Let O be the origin and let = (x^T/j). 

Proposition 2.2 One has: CCM(P) = 

/ re— 1 re— 1 

1 I 2,2 1 2 2 2 1 2 I 2 

7 a ,p\ I ^ —yiUi+i -T ViVi+l ~r %iVi+l — %i + lVi, ~ X i+lVi + + ^i^i+l " 

^ ' \i=0 i=0 
're— 1 re— 1 

\ ^ / 2 _|_2 _2 _ 2 \ \ ^ _/2 _i_2 _2 _2\ 

2_ u Vi\ X i-l * Vi-1 X i+1 Vi+l)} / j x i\ X i-l ~T~ Vi-l x i+l Vi+l, 



x i x i+l 



AA(P) 



i=0 i=0 



Proof. The perpendicular bisector to OVi has equation 

•Ei ( x i\ , Vi 

Intersecting with the perpendicular bisector to OV^+i, we find that the coor- 
dinates of the circumcenter of AOV^Vi+i are 

Xi+iVi — xjyi+i — UiUi+i + UiUi+i xjxi+i — Xixj +1 + Xi + \yf — Xiyf +l 

x — — , y = — 

2(xiy i+1 - x i+1 yi) 2(xiy i+1 - x i+1 yi) 



The denominators are four times the areas: 
2A(AOV i V i+1 ) = det(0\ OV i+1 ) = det ( %i Xl+1 ) = x iVi+1 - x i+lVi . 

V Vi Vi+1 J 

This implies the result. □ 

Remark 2.3 Rotating the Circumcenter of Mass by 90°, we get 

^ n— 1 

J^MW^-V^). (2) 



AA(P) 



i=0 



This formula is given in [2] as an integral of the recutting of polygons. 

It is also interesting to compare with the formula for the center of mass 
of a homogeneous polygonal lamina; this formula has a similar structure (the 
ratio of a cubic polynomial and the area): CM(P) = 



' n— 1 n— 1 



^ ' \i=Q 1=0 



We now use Proposition 2J2 to give a different proof of the independence of 
CCM(P) of the choice of the origin O. We wish to show that parallel trans- 
lation of the point O through some vector results in translating CCM(P) by 
the same vector. 

Let the translation vector be (£,77), that is, 

Xi = Xi + f , yi = yi + 77. 



Let us compute the formula of Proposition 2.2 in terms of The area 

A(P) is translation invariant. Let us consider the first coordinate, multiplied 
by AA(P). As a polynomial in £,77, it is: 

A + Bi + Cr] + Dfr + E V 2 , 

where A,B,C,D,E are polynomials in Xi, jji- The term A is the same as in 
Proposition 2.3; as to the rest, one easily computes: 

B = 2 ( x i-i ~ x i+i)Vh C = ^K 2 _i - x 2 +1 + y 2 _ x - y 2 i+l + 2y i (j/ i _ 1 - y i+1 )], 
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where all sums are cyclic. We see that 



B = 4A(P), C = D = E = 0, 

the latter equalities due to "telescoping" of the cyclic sums. Therefore, the 
first coordinate of CCM(P) is shifted by £. A similar computation shows 
that the second coordinate is shifted by rj. 

Next, we comment on the issue of degenerate triangles in the definition 
of the Circumcenter of Mass. A degenerate triangle is "dangerous" if its 
circumcenter is at infinity, that is, if the circumradius is infinite. This radius 
is given by the formula 

n h r 

R 



2 sin a 2 sin (3 2 sin 7 

where a,b,c are the side lengths of a triangle and a, /3, 7 are the respective 
angles. This makes it possible to distinguish between dangerous and safe 
degenerations, see Figure [3] a degeneration is dangerous if an angle tends to 
zero or 7r while the opposite side does not tend to zero. 



Figure 3: The triangle on the left is dangerous and the one on the right is 
safe 



Remark 2.4 Lemma 2.1 follows from the independence of CCM(P) on the 
choice of the point O. Indeed, choosing the point O close to vertex A on 
the bisector of the angle DAB, see Figure [4j one obtains, in the limit, two 
safe degenerate triangles OAD and OAB. With this choice of O, one has 
CCM(ABCD) = CCM(ABC © CD A), and likewise for the the other pair 
of triangles. 



3 Properties of the Circumcenter of Mass. 
Archimedes Lemma 

The center of mass satisfies the Archimedes Lemma which states the follow- 
ing: if an object is divided into two smaller objects, the center of mass of 
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Figure 4: An alternative proof of Lemma 2.1 



the compound object lies on the line segment joining the centers of mass of 
the two smaller objects; see, e.g., [3]. It turns out, the Circumcenter of Mass 
satisfies the Archimedes Lemma as well. 

Theorem 1 Let P = VqVi ■ ■ ■ V n -\ be a polygon, and let V0X1X2 . . . X m Vk be 
a polygonal line from vertex V to vertex V^. Let Q = VqV\ . . . VkX m X m _i . . . X\ 
and R = V XiX 2 . . . X m V k Vk+i ■ ■ ■ V n -\ be the two closed polygons resulting 
from cutting P along the polygonal line. Then 

CCM(P) =CCM{Q®R), 

see Figure^ 



Proof. The idea of the proof is very simple: the contribution of the triangles 
whose sides constitute the cut cancels. More specifically, fix a generic point 
O. Then 

n-l 



CCM(P) = _L_^A(AOV r i V i+1 ) CCM(AOViV i+1 ). 

Also, A(Q) CCM(Q) = 

fc-i \ 

J2M^OViV t+1 ) CCM{AOViV l+l ) +A(AOV k X m ) CCM(AOV k X m ) 

. i=0 / 
fm-l \ 

^A((AOX,X m )) CCM(AOXiX i+1 ) +A(AOX l V ) CCM(AOX 1 V ), 



i=l 
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Vk 

Figure 5: Archimedes Lemma 

and A(R) CCM(R) = 

/ m— 1 

A(AOV X 1 ) CCM(AOV X 1 ) + I ^ A{AOXiX i+1 ) CCM(AOX % X i+l ) 



i=i 

'n-l 



+A{AOX m V k ) CCM{AOX m V k ) + ^A(AOVW+i) CCM{AOV % V i+1 ) 



Adding and canceling terms, the desired equality follows. □ 

Remark 3.1 One subtelty of Archimedes Lemma has to do with degenerate 
triangulations. At first sight, a counterexample to Theorem [T] is a right 
isosceles triangle P, which is bisected into two smaller right isosceles triangles 
Q and R along its axis of symmetry. In this case CCM(P) is the midpoint 
of the base of P, which is not on the line joining CCM(Q) and CCM(R), the 
midpoints of its other two sides. 

The issue is that there is a hidden degenerate triangle. Label the isosceles 
triangle ABC with B having a right angle and let D be the foot of the 
perpendicular from B to side AC. If we perturb D away from line AC, we 
obtain an additional triangle ACD. This triangle has minute area but its 
circumcenter is far away. In this situation everything works out as expected. 
The problem occurs when we disregard this triangle because of its zero area. 
In other words, the CCM of ABC is the weighted average of the CCMs of 
ABD, BCD and ACD. 
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As a consequence of Theorem [TJ one can use any triangulation of P to 
define CCM(P), not only a triangulation obtained by connecting point O to 
the vertices. In particular, we have the following, expected, corollary. 



Corollary 3.2 The Circumcenter of Mass of an inscribed polygon is the cir- 
cumcenter. 

Proof. Let P be inscribed, and let O be the circumcenter. Consider a 
triangulation by diagonals. Then O is the circumcenter of each triangle 
involved, and hence the Circumcenter of Mass is O. 

Alternatively, one can use formula ([2]). For a circumscribed polygon, all 
\Vi\ are equal, hence ^ yields zero, that is, the origin O. Rotating back 90° 
is O as well. □ 

Remark 3.3 It follows from Theorem Q] that the Circumcenter of Mass of 
a quadrilateral is the intersection point of the perpendicular bisectors of 
its diagonals. These perpendicular bisectors may be parallel; this happens 
exactly when the area of the quadrilateral is equal to zero, see Figure [6j 



Figure 6: The perpendicular bisectors of the diagonals of the quadrilateral 
ABCD are parallel if and only if its signed area vanishes 

Next, we consider the case when P is equilateral, that is, all the sides 
have equal lengths. 

Theorem 2 If P is an equilateral polygon then CCM(P) = CM(P). 




A 



C 
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Proof. We present a computational proof based on Proposition 2.2 and 
Remark 2J3 It would be interesting to find a geometric argument as well. 

Using the formulas for CCM(P) and CM(P) from Proposition 2.2 and 
Remark 2J3 we find: 

2CCM (P) - 3CM (P) = ( ^fc/fo+i - ViVi+i ~ x lXi+l y t+l + x^y^), 

^^(xixf +1 — x^Xi + i — Xiyiyi+i + Xi + \yiyi + i)^j 

(3) 

Assume now that P is equilateral: (xj — Xj+i) 2 + (t/j — yi+i) 2 = 1 for all 
i. Then 

+ Vi + rf+i + Vi+x = 1 + 2(xiX i+ i + yiy i+ i), 

and hence 



-i + y 2 i- x - x 2 l+1 - y 2 l+1 = 2(xi-iXi + yi-xyi - XiX l+1 - y l y i+ i). 



Substitute this to the second formula of Proposition 2 J 2_ to obtain 

CCM(P) = ^Tps ( ^{ViXi-iXi + yi-xy\ - yiXiX i+1 - y 2 y l+ i), 
^(-x^x 2 - x^i-iyt + x 2 x i+ i + Xiyiy i+ i)) 

= 2A{p) i^2^ XiXi + 1 + ViVi+iliVi+i - yi)i^{ x i x i+i + yiVi+i)(xi - x i+1 )) . 

Comparing with Q, we see that CCM(P) = 3CM(P) - 2CCM(P), hence 
CCM(P) = CM(P). □ 

Another case when the Circumcenter of Mass coincides with the center 
of mass is the continuous version of the former. Let ^(t) be a parameterized 
smooth curve and O be a point not on any tangent line to 7. In other words, 
7 is star-shaped with respect to O. One defines the Circumcenter of Mass: 
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where C(t) denotes the limiting e — > position of the vector from O to the 
circumcenter of the infinitesimal triangle C"j(t)'y(t + e), and dA is the area 
of this infinitesimal triangle. Denote by CM( 7 ) the center of mass of the 
homogeneous lamina bounded by 7. 

Theorem 3 One has: CCM( 7 ) = CM (7). 

Proof. It is easiest to consider the continuous limit of formulas (|2]). Let 
j(t) = (x(t),y(t)). Then the continuous limit of Vi+i — is 27' (t), and we 
obtain 



(x 2 + y 2 )x'dt, / (x 2 + y 2 )y'dt) 



h{-II yiA JI xdA )' 



where = dxAdy and the equality is due to the Stokes theorem. It remains 
to turn this vector back by 90° to obtain 

xdA, / / ydA 



A{i) 

a well known formula for the center of mass. □ 

Thus one may view the Circumcenter of Mass of a polygon as a "different" 
discretization of the center of mass of a lamina bounded by a continuous 
curve. 



4 Generalized Euler line. Other centers 

The line through the center of mass and the circumcenter of a triangle is 
called the Euler line. The points of this line are affine combinations Ct ■= 
tCM + (1 - t)CCM; thus C = CCM and G x = CM. For example, the 
orthocenter of a triangle lies on the Euler line and is given by 3CM — 2CCM, 
that is, the orthocenter is the point C 3 . 

For a fixed t, one can repeat the construction of the Circumcenter of Mass 
of a polygon to obtain the center Ct(P): triangulate P and take the sum of 
the centers C of the triangles, weighted by their areas. The resulting center 
is again independent of the triangulation, and all these centers lie on a line 
that we call the generalized Euler line of the polygon P. Since the ratios of 
distances on the Euler line of a triangle are transferred over to those on the 
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generalized Euler line, we conclude that C t (P) = tCM(P) + (l-t)CCM(P) 
for every polygon P. 

We note that, for the case of quadrilaterals, the generalized Euler line 
was discussed in [13]. One remark we may add to the discussion is the next 
result, which follows easily from our approach. 

Proposition 4.1 Leta,b,c,d be the distances from C t (AP> CD) to C t (BCD), 
C t {CDA), Ct(DAB) andCt(ABC), respectively. Then 

aA(BCD) = cA(DAB),bA(CDA) = dA(ABC), 

see Figure^ 




Figure 7: Proposition 4.1 



Next, we shall show that the above described centers C t (P) are the only 
ones satisfying natural assumptions. Namely, assume that one assigns a 
"center" to every polygon so that the center depends analytically on the 
polygon, commutes with dilations, and satisfies the Archimedes Lemma. 

Theorem 4 A center satisfying the above assumptions is Ct(P) for some t. 
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Proof. Consider an isosceles triangle with base of length 2 and base angles 
a. By symmetry, the center is on the axis of symmetry. Let it be at height 
f(a) above the base. For example, if the center is CM then f(a) = (tan a)/3, 
and if the center is CCM then f(a) = — cot 2a. From the symmetry of the 
equilateral triangle, f(ir/S) = l/v3. 

Any triangle can be triangulated into three isosceles triangles, hence the 
function f\ot) determines the center uniquely. We shall prove that / satisfies 
a certain functional equation, solve it, and deduce the desired result. 

Consider a kite made of two isosceles triangles with angles a and j3. By 
the Archimedes Lemma, the center is on the axis of symmetry at (signed) 
height 

f(a) tana — f((3) tan/3 
tan a + tan (3 

see Figure [8j 




Figure 8: Triangle and kite 



Consider a triangle and its center, see Figure [9j One can reflect the 
triangle in either of its sides to obtain a kite. In this way, we can determine 
the signed distance from the center to each of the altitudes. For example, 
reflecting in the side AC, we use the (scaled) formula Q to find that the 
distance to the altitude from vertex B equals 

f(a) tan a — f(j) tan 7 
tan a + tan 7 

where hs — \BQ\ is the length of this altitude. 

Let cIa, cIb, dc be the signed distances from a point to the altitudes AP, BQ, 
CR, respectively. The sign is determined by orienting the altitudes from the 
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Figure 9: A center of a triangle 
vertices of the triangle. We claim that 



dg d c 
Ha h B h c 



0. 



(5) 



Indeed, the left hand side of ^ is a linear function on the plane, and it suffices 
to show that it vanishes in the vertices of the triangle. Consider vertex A. 
The triangles ABQ and ACR are similar, hence \AQ\/\BQ\ = \AR\/\CR\. 
Taking the orientation into account, this proves (|5| for point A, and likewise 
for B and C. 

Thus ^ implies the following functional equation on /: 
f(a) tana — /(/3) tan/3 f(/3) tan/3 — f(j) tan 7 f(j) tan 7 — /(a) tana 



tan a + tan /3 



tan (3 + tan 7 



tan 7 + tan a 



(6) 



for all triples of angles satisfying a + /3 + 7 = f . One can check that (]6j) 
holds for /(a) = tana and for /(a) = cot 2a. 

To solve equation (|6j), we set x = tana, y = tan/3, z = tan 7 and 
g(x) = x/(arctanx). Then (|6| becomes 



g(g) - 0(2/) + g(y) -g(z) + g(z) - gCg) = , 

x + $/ 2/ + £ ,2 + X 

subject to arctanx + arctany + arctanz = |. We have 

7T 1 — xy 

z = tan( arctanx + arctany) = . 

2 x + y 

We set y = in equation to obtain 

g(s)-g(0) | 9(0) -9(1) | g(|)-g(s) 



(7) 



(8) 



x 



+ X 
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or 

3 (x)-x 4 9 (-)+9(0)(x 4 -l) = 0. 

x 

Expand g(x) in a power series: g(x) = a + aix + a 2 x 2 + .... Then, equating 
terms of the same degrees, we find that a n = for n > 5, and = 0, 03 = a%. 

So g(x) is of the form ao + a\X + CL2X 2 + a\x 3 . Since equation (|7p is linear 
in g, it suffices to consider each monomial individually. The 0th, 1st and 
2nd-degree monomials satisfy the equation. For the 3rd-degree, we have 

fx 3 — y 3 y 3 — z 3 z 3 — x 3 \ 

at 1 1 = 0. 

\ x + y y + z x + z J 

If a\ 7^ 0, it is necessary for the expression in the brackets to vanish 
identically after substituting for z via equation ([8j. However, it is easy to 
see that it does not by, e.g., substituting x — 1, y — 2 and z = —1/3. It 
follows that ai = 0, so that g(x) = ao + d2X 2 . 

Therefore f(a) = c\ cota + c 2 tana. Since f(ir/3) = l/\/3, 

^ + c 2 V3 = ^ c 1 = l-3c 2 . 

Hence 

f(a) = (1 — 3c) cot a + ctanct = t ( ^ & — ] + (1 — t)(— cot 2a) 



with t = 3 — 6c. Thus the solution set is the Euler line, as claimed. □ 

We now show that the generalized Euler line is highly sensitive to sym- 
metries. As such, one might consider utilizing the Euler line to detect sym- 
metries, or as a certificate for their non-existence. The symmetries involved 
are of different type, some being Euclidean and some of another sort, such 
as equilateralness. 

As we saw in the previous section from Theorem |2j if P is an equilateral 
polygon, then the Euler line of P degenerates to a point. 

Theorem 5 Let P be a polygon of nonzero area so that the Euler line E is 
defined. 

1. If P has a line of reflection symmetry L, then E is either L or a point 
on L. 



17 



2. If P has a center C of rotational symmetry and none of the extensions 
of sides of P pass through the center, then E = C. 

3. Assume that the sides Si, S 2 ■ ■ ■ , S n of P satisfy \Si\ = \S 2 \ = ■■■ = 
\S n -i\. Then E is orthogonal to side S n . 

Proof. (1). First, we note that P cannot have sides which are subsets of 
L because then at least one such side would have three edges incident to a 
vertex. Thus we may place O on L and triangulate. Since every triangle 
on one side of L has a mirror image counterpart on the other side of L, the 
Circumcenter of Mass lies on L and so does the center of mass. 

(2) We let O — C and triangulate from O. To each triangle in the triangu- 
lation, there is a class of triangles obtained from it via rotation. The vector 
sum of the centers of mass of the triangles of this class is equal to O. 

(3) Let Vi, V2, . . . , V n be the vertices of P. We reflect P in side S n . Denote 
the image of a vertex under reflection by a prime. Consider the polygon 
Q = V1V2 • • ■ V n y^_i • • • V 2 - This polygon has twice the area of P and since 
it is equilateral, its Euler line is a point (the center of mass). Since Q has S n 
as an axis of symmetry, its center of mass lies on S n (or its extension). By 
Archimedes' Lemma (Theorem [T]), the Euler line of Q is the weighted sum 
of the Euler lines of P and its reflection. Since the area of the reflection is 
the same as that of P, the weights are the same. Let e and e' denote the two 
Euler lines. These are symmetric about S n , so the Euler line of Q must be 
the midpoint of every pair of mirror points of e and e'. This is only possible 
if e = e' and are orthogonal to S n . □ 

Remark 4.2 The generalized Euler line that we have defined is as a prop- 
erty of the polygon, rather than the set of its vertices. For example, if a 
different choice of edges is given to a quadrilateral, e.g., rather than consid- 
ering ABCD, we consider ABDC, its center of mass and Circumcenter of 
Mass will be different, and consequently its Euler line also. 

In the case of a quadrilateral, a definition of the Euler line which is a 
property of the set of vertices is given in |12j . following the work [16] . The 
authors of [TU] construct a point which is a "replacement" to the circumcenter 
of a quadrilateral when the quadrilateral is no longer cyclic. This point, called 
the "isoptic point", is a center of the quadrilateral which depends only on 
the set of vertices, in the sense that it is independent of the choice of edges 
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for the four vertices of the quadrilateral, e.g., it is the same for ABCD and 
ABDC . The Euler line is then defined as the line connecting the isoptic 
point and the centroid of the vertices of the quadrilateral. In [18] , the second 
author extends the definition of the isoptic point to (n + 2)-polytopes in M. n , 
thus constructing an Euler line which depends only on the set of vertices of 
the (n + 2)-tope. 



5 In higher dimensions 



In this section, we extend the construction of the Circumcenter of Mass to 
simplicial polyhedra in W 1 . The construction is the same: given an oriented 
simplicial polyhedron P, choose a point O that does not belong to any of the 
hyperplanes of its facet, and triangulate P by the simplices with vertex O 
and the bases the facets of P. Let C{ be the circumcenters of these simplices, 
and define the circumcenter 

Vi 



CCM(P) = J2V(P) Q 



where V, denotes the signed volume of ith simplex, V(P) is the volume of P, 
and the sum is taken over all facets of P. 

Similarly to the 2-dimensional case, we have the following result. 

Theorem 6 The Circumcenter of Mass is well defined, that is, does not 
depend on the choice of point O. 

As a preparation to the proof, we state two lemmas. The first one is well 
known and we do not prove it. Let S = (Vq, V±, . . . , V n ) be a simplex in M. n . 

Lemma 5.1 The signed volume Vol (S) is 

detMQg) 
n\ 

where M is an n + 1 by n + 1 matrix 

V Q Vi ... V n 
1 1 ... 1 

The second lemma is lesser known so, for completeness, we provide a 
proof. Let S = (Vq, Vi, . . . , V n ) be a simplex and C(S) its circumcenter. 
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Lemma 5.2 The ith coordinate of C(S) is 



detMj(g) 
2detM(S)' 

where the matrix Mi is obtained from the matrix M by replacing its ith row 
by the row (|F | 2 , |^i| 2 , • • • , |K| 2 )- 



Proof. Let V$ = (xn, . . . , x in ), i — 0, . . . , n. Consider the equation 

/y»2 I m m I rp 2 /y»2 I ^ I ^v»2 ,y»2 I ^ ^ ^ I ,yi 2 ™ 2 I ... I /y»2 



det 



X'2 
r j. 



£01 

X 2 
1 



ari2 

1 



■Erin 
1 



in variables x%, . . . , x ra . This is an equation of a sphere, and this sphere passes 
through points V , . . . ,V n . Thus this is the equation of the circumsphere of 
the simplex S. 

Expand the determinant in the first column: 



(x\ + • • • + x 2 n ) det M — x\ det M x - x 2 det M 2 



x n detM n + C = 0, 



or 



Xi 



det M x 
2detM 



H +[x 



detM n 
2detM 



Ci, 



where C, C\ are constants depending on the coordinates of the vertices of the 
simplex. This implies the statement of the lemma. □ 



Proof of Theorem [6} Consider the quadratic map Q : W 1 — > R given by 
the formula Q{V) = \V\ 2 . The graph of Q is a paraboloid. For V G M n , let 
V e W n+1 be the vector (V, Q{V)). If S is a simplex in R n , let 5 be its lift 
in E n+1 . 

Consider the linear map 7Tj : ]R n+1 — > W 1 given by the formula 

(Xi, . . . , X n +i) I y \Xij ■ ■ ■ , 2Tj_i, X n j t \^ 2Tj_|_i, . . . , Xn). 

In terms of these maps, the ith coordinate of C(S) is Vol (7Tj(5'))/2Vol (S). 
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Let F be a facet of P. Denote by Sp the cone over F with the vertex O. 
The statement that we need to prove is that, for every i, 

Vol (9) 

FcdP 

is independent of O, where the sum is taken over the facets of P. Without 
loss of generality, consider the case i = 1. 

Let Vi, . . . , V n be the vertices of a facet F. The respective summand is 

\0\ 2 IVil 2 ... IKP 



det 



... V n 

1 1 ... 1 



(10) 



where O = (o 2 , . . . , o n ) and likewise for Vi, . . . , V n . 

The coefficients of \0\ 2 , o 2 , o 3 , . . . , o n in (10) are, up to the sign, 

Vol (p 1 vr 1 (F)),Vol (p27n(F)),...,Vol (^vn(F)), 

where the projection pi : R n — )■ IR n_1 is given by forgetting ithjx>ordinate. 

It remains to notice that the total signed volume of pi%i{dP) is equal to 
zero since OP is a cycle. Therefore the coefficients of \0\ 2 , o 2 , 03, . . . , o n in 
(P) vanish, and this sum does not depend on O, as needed. □ 



Similarly to Proposition 2.2, one has an explicit formula for CCM(P) 
that follows from Lemmas 5^ and 52 and the proof of Theorem [6] 

Let F = (Vi,...,V n ) be a face of a simplicial polyhedron P C M. n . 
Let A(F) be the matrix whose columns are the vectors Vi, . . . , V n , and let 
Ai(F) be obtained from the matrix A(F) by replacing its ith row by the row 
(| Vi\ 2 , . . . , |V n | 2 ). Choosing O as the origin yields the following formula for 
the ith component of the Circumcenter of Mass. 

Proposition 5.3 One has: 

CCM(P)i = — — - V det AAF). 
v ' 2(n\)V(P) ^ 

v ' v 1 FcdP 

This rational function of the coordinates of the vertices of the polyhedron 
P has homogeneous degree one: it is a ratio of polynomials of degrees n + 1 
and n. 



As an application of this formula, let us prove an analog of Corollary 3.2 



Corollary 5.4 The Circumcenter of Mass of an inscribed polyhedron is its 
circumcenter. 
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Proof. Choose the origin O at the circumcenter and, without loss of gen- 
erality, assume that the circumscribing sphere is unit. We wish to prove that 
CCM(P) = O; to be concrete, consider the nth component of CCM(P), as 
given in Proposition |5.3| 

The matrices A n (F) have nth row consisting of Is, that is, these are the 



matrices from Lemma 5.1 It follows that, up to a factor, CCM(P) n is the 
signed volume of the projection of dP to M n_1 along the last coordinate 
direction. This volume vanishes, as in the proof of Theorem [6j and we are 
done. □ 

A multi-dimensional version of Theorem [T] holds as well, and for the 
same reason: the contribution of the simplices whose bases belong to the cut 
cancel out. Thus, with the self-explanatory notation, one has the Archimedes 
Lemma. 

Theorem 7 If a simplicial polyhedron P is decomposed into simplicial poly- 
hedra Q and R then CCM(P) = CCM(Q © R). 

It follows that one can find the Circumcenter of Mass using any triangu- 
lation of P. 

Similarly to the 2-dimensional case, one can take an affine combination of 
the centroid and the circumcenter of a simplex to define the center Ct(P); as t 
varies, these centers lie on a line that we call the Euler line of the polyhedron. 

The Euler line of a simplex is well studied: this is the line through the 
centroid and the circumcenter, see IB] and the references in the latter 
article. 

The Euler line of a simplex contains its Monge point. This is the intersec- 
tion point of the hyperplanes through the centroids of the n — 2-dimensional 
faces of a simplex in IR n , perpendicular to the opposite 1-dimensional edges 
(in dimension two, the Monge point is the orthocenter). In our notation, the 
Monge point is C( n+ i)/( n _i). Thus we obtain a definition of the Monge point 
of a simplicial polyhedron in W l . 



6 On the sphere 

In this section, we consider a version of the theory on the sphere (and, to 
some extent, in the hyperbolic space). Let us start with S 2 . 
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A circle on the unit sphere in M. 3 is its section by an affine plane, that is, 
a plane not necessarily through the origin which is the center of the sphere. 
To an oriented circle we assign its center; the choice of the two antipodal 
centers is made using the right-hand rule. Given a spherical triangle ABC, its 
circumcircle is oriented by the cyclic order of the vertices, and this determines 
the circumcenter of the oriented triangle. 

We need the notion of the center of mass of a collection of points on the 
sphere, see [7J. Let Vi G S 2 , i = 1, . . . , n, be points and rrii > their masses. 
The center of mass of this system of masses is the point 

\Y, m i V i\' 

taken with the mass | '^2 i m i Vi\. This notion of the center of mass satisfies 
natural axioms and is uniquely characterized by them [7J. 

The notion of the center of mass naturally extends to spherical lamina. 
If U is a spherical domain then its center of mass is the vector 

v dA, 



u 

normalized to be unit, where v G S 2 is the position vector of a point and dA 
is the standard area form on the sphere. 

In particular, let ABC be a spherical triangle. Then the center of mass 
of the triangular lamina is the vector 

AxB Mi3 vBxC rf (^) ,CxA d{C ' A) (11) 
AXB sin d(A,B) +BX C sind(B,C) +C A sin d(C, A)' {U) 

normalized to be unit. Here Ax B is the cross-product in IR 3 , and d(A, B) is 
the spherical distance between points, that is, the side length of the triangle. 



The mass of the triangular lamina is the norm of the vector (11). Formula 



([TT]) can be found in U U 03] • 

With these preparations, we can define the circumcenter of mass of a 
spherical polygon and study its properties, similarly to the Euclidean case. 

Let P = (Vi, . . . , V n ) be a spherical polygon and let W G S 2 be a point, 
distinct from the vertices of P. Consider the triangles WViVi + \ where the 
index % is cyclic. Let Oj be the circumcenter of ith triangle, taken with the 
mass rrii, equal to the area of the plane triangle WViVi + \. Define CCM(P) 
to be the center of mass of the collections of points (Oj,raj). 

Theorem 8 The point CCM(P) does not depend on the choice ofW. 
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Proof. Let ABC be a spherical triangle. Its circumcenter is the intersection 
of the sphere with the ray, perpendicular to the plane of the triangle and 
oriented according to the orientation of the plane given by the cyclic order 
of the vertices. Hence the circumcenter is the vector (A — B) x (A — C), 
normalized to be unit, that is, 

AxB+BxC+CxA 



\AxB + BxC + CxA\ 



Note that the denominator is twice the area of the plane triangle ABC. 

Using this formula for the triangles WViVi + \ and taking the weighted 
sum, yields 

CCM(P) = (12) 

I V i X V i+l I 

since x V i + V i+i x W) = 0. Thus CCM(P) is independent of W. □ 

The Archimedes Lemma holds, the same way as in the Euclidean case. 
We also have an analog of Theorem [2j 

Theorem 9 If P is an equilateral spherical polygon then its circumcenter of 
mass coincides with the center of mass of the lamina bounded by P , that is, 
CCM(P) = CM(P). 

Proof. Let £ be the side length of P. Choose a point W and triangulate 
the polygon as above. Let (Qi, rrii) be the center of mass of the triangular 
lamina WViVi + i. Then CM(P) is the center of mass of the system of points 



(Qi,rrii). By formula (11), this is 



V^r, rr * V^tt, rr d{W,Vi) „, d(W, V i+1 ) 

^— ' sm£ sina(vv,Vi) smd{W,Vi +1 ) 

as always, normalized to a unit vector. The last two sums cancel each other 



and we obtain CCM(P), as given by formula (12). □ 



An analog of Theorem [3] holds as well. Let 7(t) be an oriented simple 



closed spherical curve. The continuous limit of the sum (12) is the vector- 
valued integral 

7 x 7' dt 
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(note that it does not depend on the parameterization). Denote this integral 
by CCM(7), and let CM(^) be the center of mass of the homogeneous lamina 
bounded by 7. 

Theorem 10 One has: CCM( 7 ) = CM( 7 ). 

Proof. Let U be the domain bounded by 7. We claim that 



where, as before, v is the position vector of a point. 

Let £ be a test vector. Then what we need to establish is the equality 



Since v is the unit normal vector to the sphere, the left hand side is the flux 
of the constant vector field £ through U . The right hand side equals 



The last equality is easily verified by a direct computation, writing 7 = 
(x,y,z), £ = (a, b, c), where x,y,z are the Cartesian coordinates in IR 3 and 
a, b, c are constants. □ 

The case of higher-dimensional sphere S n C IR n+1 is similar. The cross- 
product is replaced by the skew-linear operation 






The result will follow from the Stokes formula once we show that 



curl( 7 x f ) = f . 



(Vi, 



. , K) H> Vi x 



X Vr 



n 



defined by the equality 



det(Fi,. 



V n ,0 = (V 1 x 



x K) • £ 



for every test vector £. 



25 



Given a simplicial spherical polyhedron P, one triangulates it by choosing 
a point W and takes the weighted sum of the circumcenters of the respective 
simplices. This results in the vector 

J^Vi x • • - x 14, (13) 

normalized to be unit, where the sum is taken over the facets (Vi, . . . , V n ) of 
the polyhedron P. 

Finally, we make a brief comment on the hyperbolic case. In this case, the 
Euclidean space is replaced with the Minkowski space, the pseudo-Euclidean 
space of signature (1, n), and the sphere S n by the pseudo-sphere, the hyper- 
boloid model of the hyperbolic space H n . One can still define cross-product, 



and one can use formula (13) for the circumcenter of mass. The center of 
mass is also defined in the hyperbolic case, see [?]. 

However, the geometrical interpretation is not as straightforward as in 
the spherical case. To fix ideas, consider the case n = 2, the hyperbolic 
plane. 



One has three possibilities for the vector (13): it may be space-like, time- 
like, or null. In the first case, the line spanned by this vector intersects the 
upper sheet of the hyperboloid z 2 — x 2 — y 2 = 1, and we obtain a point 
of the hyperbolic plane. In the last case, we obtain a point on the circle 
at infinity, and in the second point on the hyperboloid of one sheet 

z 2 — x 2 — y 2 = — 1 (outside of our model of H 2 ). 

This issue is already present in the case of a triangle. The intersection 
of the upper sheet of the hyperboloid z 2 — x 2 — y 2 = 1 with an affine plane 
may be a closed curve, disjoint from the null cone z 2 = x 2 + y 2 ; it may be a 
closed curve tangent to the null cone; or it may be an open curve intersecting 
the null cone. A triangle is circumscribed by a curve of constant curvature, 
the intersection of the hyperboloid z 2 — x 2 — y 2 = 1 with the plane of the 
triangle. If the curvature is greater than one, this curve is a circle; if the 
curvature equals one, it is a horocycle; and if the curvature is less than one, 
the curve is an equidistant curve. These three cases correspond to the three 
relative positions of the plane and the null cone. 

All three types of curves are represented by circles in the Poincare disc 
model. In the first case, the circle lies inside the Poincare disc, in the second 
it is tangent to the boundary, and in the third it intersects the boundary. 
Only in the first case the center of the circumcenter of a triangle is a point 
of the hyperbolic plane. 
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Remark 6.1 The results on the circumcenter of mass in the Euclidean case 
can be obtained in the limit R — >■ oo from the spherical case, where R is the 
radius of the sphere. 

Remark 6.2 A proof that the Bicycle (Darboux) Transformation [T7] pre- 
serves the CCM of triangles and quadrilaterals in spherical and hyperbolic 
geometry is entirely analogous to the one in Euclidean space. This fact and 
conservation of area under the Transformation allows to fully describe the 
dynamics of these polygons under the Bicycle Transformation. 
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